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Abstract 

The response to a magnetic flux is considered of the vacuum state of charged Dirac 
fermions interacting with a domain wall made of a neutral spinless field in (3+1) di- 
mensions with the fermion mass having a phase variation across the wall. It is pointed 
out that due to simple C parity arguments the spontaneous magnetization for this 
system is necessarily zero, thus invalidating some claims to the contrary in the litera- 
ture. The cancellation of the spontaneous magnetization is explicitly demonstrated in 
a particular class of models. The same calculation produces a general formula for the 
electric charge density induced by the magnetic flux - an effect previously discussed 
in the literature for axionic domain walls. The distribution of the induced charge is 
calculated in speciflc models. 



1 Introduction 



The possibility that domain walls that could have existed in the early universe could also 
be related to generation of a primordial magnetic field correlated at large distances [|I| has 
been recently discussed in the literature [2-5]. The discussed models are based on the idea 
that fermions, coupled to the field forming the wall, develop a spontaneous magnetization 
perpendicular to the wall. Although it is a simple exercise in general physics to show that 
a uniform magnetization of an infinite domain wall does not produce magnetic field @], the 
phenomenon of magnetization of the wall is interesting on its own. 

The claims to a non-zero magnetic moment of the ground state of a fermion field coupled 
to the wall in (3+1) dimensions are inferred from the behavior in a (2+1) dimensional QED 
of Dirac fermions with a definite sign of the mass term m. Namely, in certain calculations 0, 
1^, ^ of the energy of the ground state of the fermion field in an external magnetic field B, it 
is claimed that the total energy contains a linear in magnetic field term proportional to mB, 
which corresponds to a spontaneous magnetization proportional to the mass parameter m 
(including the sign). For a Dirac fermion field coupled to a domain wall in (3+1) dimensions 
the quantization of the motion perpendicular to the wall splits the fermion system into an 
infinite set of modes, each corresponding to a (2+1) dimensional QED with its own parameter 
m. If the phase of the fermion mass varies across the wall, the set of positive values of m 
differs from that of the negative values of m. Therefore in this picture it might at least be not 
obvious that the overall magnetization cancels after summation over modes corresponding 
to positive and negative values of m. 

It is nevertheless quite easy to argue that the cancellation necessarily takes place and the 
magnetization of the fermion field ground state at the wall in (3+1) dimensions is strictly 
zero. Indeed, the Lagrangian density for the fermions with the phase of the mass /z depending 
on the coordinate z perpendicular to the wall can be generally written as 

L = ^[i{d^-t A^) 7" - fii{z) - ifi2{z) 7'] ^ (1) 

with 111 and /X2 being respectively the real and imaginary parts of ii{z), and A standing for 
the electromagnetic field potential and absorbing the charge e in the normalization of the 
field. The variation of fi2 breaks the P and CP parities. However the C parity is manifestly 
conserved with both /ii and /i2 being C-even. Therefore after the fermions are integrated out 
the energy of the system as a function of A in the C-even background of fi{z) cannot contain 
odd powers of the C-odd field A. In particular the energy cannot be linear in the magnetic 
field B. The same argument holds also for the situation where the fermions are assumed to 
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^^ = — :t—-2 — E-B. (2) 



possess anomalous magnetic and/or electric dipole moment since both such interactions also 
conserve the C parity. 

This simple general C parity argument is clearly sufficient for excluding the possibility of 
spontaneous magnetization of the fermion vacuum in a domain wall background. However 
due to the existence in the literature of claims to the contrary ^ it is quite instructive to 
demonstrate explicitly the vanishing of the linear in the magnetic field Bz = B term in the 
energy (as well as of all the odd terms) at least in a specific model. Moreover the C parity 
argument certainly allows dependence of the energy on even powers of the electromagnetic 
field, which gives rise to the most interesting phenomenon of appearance of electric charge 
density once a magnetic flux is applied across the wall. This phenomenon is related to the 
well known coupling of a pseudoscalar field to the electromagnetic invariant E ■ B. For 
slowly varying with z mass term, one can approximate around a given point Zq as 
lji{z) = fi^zo) + Sfi{z), and treating the varying part 6fi{z) as small perturbation, find the 
term proportional to E ■ B in the energy density w from the well known triangle graph as 

1 /i2 Sfii - Hi 6 fl2 

47r2 /if + ^2 

The charge density is then found from variation with respect to the potential Aq. For the 
magnetic field in the z direction, Bz = B, one finds 

V-O^ hVW = -Tl- = T^:rarctan . 3 

The total charge is then given by the total flux of the magnetic field through the wall, 
F = J dxdy Bz and the difference of the phases of the mass term fj,i + ifj,2 at two infinities 
in z, A$ = arctan(/i2//ii)|2^+oo - arctan(/i2//ii)|2^-oo, 

Q = J dxdydz{^{z)^°^{z)) = ^-^. (4) 

The induced charge is a direct analog of the (generally fractional) fermionic charge of 

a kink in (1+1) dimensional models [llO| , |lT|. This phenomenon in (3+1) dimensions was 

considered [|T2| in connection with a magnetic monopole (dyon) traversing an axionic domain 

wall, in which process the net change (reversal) of the magnetic flux across the wall results 

in charge exchange between the dyon and the wallQ 

^The paper contains a reference to an unpublished communication with H. Georgi and J. Polchinski, 
who apparently had also interpreted the effect of the triangle graph in terms of an induced charge in a 
magnetic field. However, as they both kindly corresponded to me, they had never pursued this issue beyond 
an unwritten remark. The monopole - axion wall charge exchange was also later discussed by I. Kogan [0. 
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In the specific model, considered in the present paper, the calculation of the induced 
charge density automatically comes along with the calculation of the (eventually vanishing) 
spontaneous magnetization. It will be shown that, as expected on general grounds, the 
relation between the total induced charge and the total magnetic flux does not depend 
on the specific shape of the dependence of the phase $ on z. However the distribution of 
the induced charge density in z does depend on the specific rate of variation of the mass 
parameter ^{z), and generally differs from that given by the equation (§) that is justified 
in the limit where the rate of variation can be considered as slow. One can notice in this 
connection that the charge distribution may be of a greater physical relevance than the total 
charge, since barring the existence of monopoles it is physically impossible to produce a net 
magnetic flux through an infinite or closed wall. Therefore the total charge has to be zero, 
while the distribution of the density of charge can be nontrivial. 

The further material in this paper is organized as follows. In Sec. 2 the considered class 
of models is described as well as some properties of the relevant operators corresponding 
to the motion in the x — y plane in magnetic field and to motion in the z direction in the 
domain wall background are discussed. In Sec. 3 the dependence on B of the energy of the 
ground state of the fermion field is calculated and the vanishing of all odd terms in the 
expansion in B is demonstrated including the vanishing of the spontaneous magnetization. 
In Sec. 4 a general expression for the induced charge density is presented, and Sec. 5 contains 
the calculation of this density in the limit of slowly varying phase as well as a discussion of 
the topological nature of the total induced charge. An explicit calculation of the distribution 
of the induced charge in two particular sample situations is considered in Sec. 6. Finally, a 
general discussion and a summary of results are presented in Sec. 7. 



2 The model and the relevant operators 

The simplifying assumption in the class of models to be considered here is that the real part 
fii of the mass term is fixed and nonzero, while the imaginary part /i2 depends on z. To the 
best of the author's knowledge, a model of such type with fi2 being proportional to a spinless 
field 0, fi2 = g 4>, with the field varying as (j){z) = v tanh.{m^ z/2) across the domain wall, 
was first suggested by T.D. Lee and G.C Wick [|I^] as a model of spontaneous breaking 
of the CP symmetry. The fermion spectrum and scattering states in the presence of the 
domain wall in this model were studied in Ref. |jl5|. 



In a generic case of a non-zero change of fi2 between the infinities in z and a non-zero 
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constant /ii, one can assume for definiteness without further loss of generahty that jJi is 
positive, and fi2{z) changes from a negative value at z ^ — oo to a positive one aX z +00. 
Also for definiteness it is assumed here that a uniform positive magnetic field B is applied 
in the z direction. In what follows the gauge for the electromagnetic field is fixed in a 
standard way such that the vector potential for the field B is given by = 0, Ay = B x. 
Adopting also the standard representation for the 7 matrices, the one-particle Hamiltonian, 
corresponding to the Lagrangian ([^) takes the following form 



H= , (5) 










iP^ 


-iR 
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iP 






-iP 


-iR 









V 











J 



and has the diagonal form: 

= diag (/i? + ptp + jiji^^ ^2 ^ ppt + ptp^ ^2 ^ ppt + ppt^ ^2 ^ ptp ^ ptp) (Q) 

These formulas make use of the following notation for the operators describing respectively 
the Landau quantization of the motion in the x — y plane and the quantization of motion 
along the z axis: 

R = d^ + Bx+Py, R^ = -d^ + Bx+py , (7) 



and 



P 



d,+fl2{z), P^ = -d, + fM2{z) . (8) 

The quantity Py is the value of the conserved momentum in the y direction (as a consequence 
of the chosen gauge condition). The energy levels, determined by the eigenvalues of R'^R and 
RR^ do not depend on py and the degeneracy number is well known to be given by B S/ {2tt), 
where 5* is the normalization area in the x — y plane. In what follows this degeneracy factor 
will be explicitly accounted for, and the value of Py set to py = in the definition of the 
operators (0). 

According to Eq.(|^) the spectrum of one-particle energies is determined by eigenvalues 
of the operators R'^R, RR\ P'^P, and PP^ . These spectra exhibit two separate structures 
found in supersymmetric quantum mechanics, well known for the operators R and i?^ and 
also recently used for the P type operators in connection with kinks in (1+1) dimensional 
models [|16| . Namely, if the issue of boundary conditions in z is ignored, one would naively 
conclude that the spectra of eigenvalues of P'^P and PP^ coincide except for an extra zero 
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eigenvalue (under the adopted sign conventions) of P^P. Indeed, let Vk, A; = 1, 2, . . ., be 
the normalized eigenfunction corresponding to the (necessarily positive) eigenvalue of the 
positive operator PP^, so that 

PP^Vk = Xlvk. (9) 
Applying the operator to both sides of this equation, one finds that the function 

Uk = P^Vk/Xk (10) 

is the normalized eigenfunction of the operator P^P with the same eigenvalue A^. Applying 
the operator P to both sides of the latter relation and using (^, one find the inverse of the 
relation (p!0|) 

Vk = Puk/Xk ■ (11) 

This construction does not work however for the zero mode uq of P^P, satisfying the equation 
P Uq = (BPS state). The explicit form of the normalizable function Uq{z) is readily found 
from the definition (|^): 

Uq{z) = const ■ exp ( — / /i2(5) dz ) . (12) 



The operators and R coincide, up to normalization, with the creation and annihilation 
operators for a harmonic oscillator, and their spectra are the text-book ones: the spectrum 
of eigenvalues of RR^ is given by 2 P n with n = 1, 2, . . ., while that of R^R is given by the 
same simple expression, however also including n = 0. 

This discussion of the properties of the operators involved in the Hamiltonian in Eq.(|^) 
is helpful in considering the spectrum of the one-particle energies. In particular, one can 
separately consider each eigenmode of the motion in the z direction as a (2+1) dimensional 
fermion system. Then the eigenmodes of the operator P^P correspond to such system with 



positive mass parameter m: m = yP'^P + fif, while the eigenmodes of PP^ correspond 
to negative m: m = —JPP'f + /if. For the former modes the negative energy spectrum 



of Landau levels is given by — wPtP + /if + 2 P n and includes n = 0, while for the latter 



ones the negative energies are given by — yPPt + /if + 2Pn excluding n = 0. Since all the 
non-zero eigenvalues of P^P and PP^ coincide, one might very naively conjecture that their 
effects in magnetization cancel and the net result is given by only one 'unpaired' zero mode 
of the operator P'^P. Such conjecture however would be false, since the summation over the 
modes is generally divergent, and one should perform a proper calculation with a proper 
regularization. Also the discretization of the continuum spectra of the operators P^P and 
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PP^ requires imposing conditions at the boundaries of large but finite bounding box in z. 
These conditions generally are not satisfied by the relations ( jlOD and (jTl]), and there arises 
a splitting of the spectra of continuum modes of the operators P'^P and PP'^ which should 
be accounted for in a proper calculation. 



3 Dependence of the fermion vacuum energy on B. 



For a full regularized calculation of the energy of the ground state of the fermion field we use 
here the standard four dimensional technique with the Pauli-Villars regularization procedure. 
The latter regularization preserves the gauge invariance and, importantly for the discussed 
problem, the Lorentz covarianceQ. This amounts to introducing the regulator fermion field 
with large but finite mass, so that the mass term for \E' can be written as M + i/i2 7^ instead 
of yUi + z/i2 7^ for the 'physical' fermion ?/^, and treat the loop with the regulator field with 
an extra minus sign^ The regularized expression for the total energy W then reads as 



W 



- fTr In [iP'a7° - /"i - ifJ'2iz) 7' 



Tr In 



iD,7"-M-2/X2(^)7']} 



Or 



M 



Tr / dm 

00 27r J 



iDa'j" - m - ifi2{z) 7' 



(13) 



where iDa'y'^ = 7V0 ~ 7 ■ (p + A) with the spatial momentum p understood as operator, 
and the trace running over the spinor indices and the spatial variables. 

The inverse of the Dirac operator in the last expression in Eq.(|13D can be readily found 
by the usual multiplication of the numerator and the denominator by iD^'j"' + m — ifi2{z) 7^. 
After taking into account the previously mentioned degeneracy in the momentum py, the 
result can be written as 



W 



.BS 
'"27 



dpi 



Or 



Tr / dm 

00 27r J 



M 



Po + m 


-iP 





Po + m 
-iR 
-zPt 



-iP^ 
-iR^ 
-po + m 




iR \ 
-iP 


-po + m J 



H\m) 



-1 



(14) 



^The vanishing of spontaneous magnetization can in fact be viewed as due to the possibihty of using a 
regularization preserving the Lorentz symmetry. 

•^Only one regulator field is indicated here for simplicity of expressions, whereas the full regularization of 
the vacuum energy requires additional regulator fields. However for the considered here effects in the energy 
one regulator is in fact sufficient. 
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where H'^{m) is the same diagonal matrix as in Eq.(|^) with /ii being replaced by m. The 
trace over the spinor variables leaves only the contribution of the diagonal factors in Eg. (JL4\) , 
which thus can be rewritten as 



W 



.B S r°° dpo 
2ti J-oo 2tt 



M 



dm < Tr 



Ml 



PoT 



+ Tr 



m 



Po 



(15) 



The first term in the curly braces is manifestly odd in po; thus the integration over po from 
— oo to oo yields zero. (The integral can be verified to be finite, thus no further regularization, 
that potentially could invalidate the symmetry argument, is needed.) The second term after 
explicitly performing the trace over the spinor variables takes the form 



W = -i- 



BS 



27T J-oo 27T 
1 



dpo f^' ^ 
/ dmm ir 



1 



+ 



1 



+ 



-pI + m? + PW + RW -pI + m2 + Ptp + R]R 
1 



pl + m? + PP^ + RW -pl + m? + PP^ + R'^R^ 



(16) 



Here in braces are grouped together the terms with the same ordering of P and and 
different order of R and P^. Each of the two expressions in the braces has the general form 



Tr 



1 



(17) 



.X + ppt X + ptp, 

where the operator X does not depend on the magnetic field. It is now easy to show that the 



expression (|T^ contains only odd powers in its expansion in B and thus that the expansion 
in B of the vacuum energy described by Eg . (p!6|) contains only even powers, as reguired 
by the C invariance. Indeed, applying the formula = exp{—i3 x) djS and using the 
described spectra of the operators PP^ and P^P, one can perform the trace over the space 
of the latter operators and find the expression (p!7| ) in the form 

dp coth (PB) Trexp (-/?X) , (18) 

which is manifestly odd in P0. 

This concludes the proof by an explicit calculation of the absence of odd powers of the 
magnetic field B in the expansion of the energy of the fermion vacuum. 

^Clearly the first term of the expansion of the expression in Eq.(p^, proportional to B~^, results in a 
divergent expression, which corresponds to the divergence of the vacuum energy at -B = 0. This is the only 
place where additional Pauli-Villars regulators are formally required. The subsequent terms of the expansion 
however give a finite difference W{B) — W{Q) that is even in B. 
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4 Fermion charge density induced by magnetic field 



We now proceed to a calculation of the electric charge induced by the external magnetic field 
B. The general formula for the charge density p{z) = {'ip(z) 'j'^ ^{z)) is obtained from the 
variational derivative of the energy in Eq.(|l3|) with respect to Aq{z) at = and reads as 

pW = -^/^^(^|Tr{7o [^D^r-^^,-^^^,{z)J']'']\z'^ . (19) 

The charge density is finite thus eliminating the need for considering the contribution of the 
Pauli-Villars regulator. The trace here runs over the spinor indices and over the space of 
the operators R and . The z dependence is thus left 'untraced', as is indicated in Eq.(|T^) 
by the diagonal matrix element of the remaining (after 'tracing out' the spinor variables 
and the motion in the x — y plane) z dependent operator in the coordinate representation. 
Performing also the trace over the z dependence, i.e. calculating the integral over z, would 
yield the total induced charge Q. Using the inverse of the Dirac operator as in Eq.(|l4D one 
can calculate the trace over the spinor indices and rewrite the latter formula for p as 



Pi 



B 
2^ 



dpo / 

271 \ 

1 



Tr 



1 



1 



^0 + /if + Ptp + R\R 
1 



Po + Pl 



ptp + RW 



,Po + Pi 



ppt + ptp pI + pI + pp^ + ppt 



(20) 



In the latter expression also the Wick rotation po ~^ "ipo is done and the part of the integrand 
that is odd in po is discarded. The trace over the space of the operators R and R^ in Eq. (^) 
is finite and trivial: 



Tr 



1 



.x + ptp 

Using this formula in Eq.(|20|) one finds 



x + ppt 



1 

X 



B 



^ 2tt J-oo 2tt 



dpo 



Po + Pi 



ptp pl + pi + PP^ 



(21) 



(22) 



One can notice the absence of higher than linear terms in the expansion of the charge density 
in powers of B, which is a consequence of the relation (^). It can be also observed that up 
to the factor B/ (27r) the expression ( [2^ ) coincides with similarly calculated fermionic charge 
in a (1+1) dimensional theory with z dependent mass term pi + i fi2{z) 7^. This relates the 



present calculations to the problem of fermion charge of kinks fl^, [T| 
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5 Charge density for slowly varying phase and the topo- 
logical nature of the total charge 

We now consider few specific cases of the dependence of /i2 on z, in which the charge density 
can be calculated from Eq.(P^ either in a closed form, or at least 'in quadratures'. Also the 
topological nature of the total charge is to be addressed in this section. 

The first case to be considered is the one where the rate of variation of the parameter /i2 
can be considered as slow. In this case in order to calculate the charge density at a point zq, 
one can approximate fi2{z) near zq by the linear dependence: /i2(^) ~ fJ'2{zo) + (z — zq) h^^zq) 
and also expand the difference of the Green's functions in Eg. (^21) to the linear order in 
/i2(2o): 

B dpo 



p{zo) = 2/ii/i^(zo) 7^ / -^(zo 

ZTT j-oo Zn \ 



2/ii/i2(Zo)^ y^^-^^-p- [p^ + Po + f^l + f^2{Zo) 

B /.i/.'2(^o) _ B d^^^^J^,,{z) 



At:"^ Hi + hI^Zq) 47r2 dz \ pi 



(23) 



2 = 20 



where the final expression matches that in Eq.(^) in the considered case of a constant /ii. 
Clearly, for this calculation to be justified the characteristic length of variation of P2{z) 
should be much longer than the 'local' Compton wavelength (/if + ii'^{zq))~^/'^ . 

Let us address now the statement that the total integral over the density, Q = /f^ p{z) dz, 
is in fact topological, i.e. it is determined only by the limiting values of fi2{z) at the infinities 
and does not depend on the specific shape of the function fi{z) at finite z. From Eq.(p2|) the 
total integral can be written ClS db trace: 




p§ + /if + PPt^ 



(24) 



In the latter expression the traces of the individual operator terms in the braces are 
infinite. For this reason one cannot calculate them separately, each in its own basis. If 
each trace were finite, one could calculate the trace of the first operator in the basis of the 
eigenfunctions Uk of the operator P^P and the trace of the second operator in the basis 
of the eigenfunctions Vk of the operator PP'I' (cf. equations - |ll])). In that case the 
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difference would arise only from the 'extra' zero mode (|T2D of the operator P^P. In view 
of the divergence of each of the traces only the trace of the whole operator in braces in 
Eq.(p^ can be calculated in an arbitrary complete basis. Also the continuum spectrum of 
the operators should be discretized by choosing a large bounding box in z: conventionally 
from —L/2 to where L is sufficiently large for the asymptotic values of ^2{z) to set 

on, and by imposing boundary conditions at z = ±L/2. Consider, for instance, the case 
where the trace in Eq. (p^ is calculated in the basis of the eigenf unctions Vk of the operator 
PP^ with an (anti)periodic condition at the boundaries of the box. Then the eigenf unctions 
Uk of the operator P^P found from Eq.(|lOD do not satisfy the same condition since fi2 takes 
different values at the infinities. This effect thus gives rise to a splitting of the eigenvalues 
in the two terms in the braces in Eq. (|2^) . The latter effect however vanishes, if fi2{z) takes 
the same value at both infinities. Also in this case the operator P'^P no longer has the zero 
mode described by Eq. ([12|) , and the spectra of the two discussed terms in Eq.(|^) coincide 
in the same basis of (anti)periodic eigenfunctions. Thus the total charge vanishes, if there is 
no net variation of fi2- /i2(— oo) = ^12(00). 

The latter observation in combination with Eq.(P5|) is in fact sufficient to conclude that 
the total charge depends only on the difference of the limiting values of fi2- Indeed, consider 
a large box, defined by the length L, and let the difference A = fi2{L/2) — L/2) be 
nonzero. Consider now an extended system, where the box is continued from z = L/2 to 
z = L, in such a way that the length of the extension L — L/2 is arbitrarily large, and the 
behavior of fi2{z) in the extension is chosen 'by hand' such that it is a slow function smoothly 
interpolating ^2 from its fixed value at z = L/2 back to fi2{—L/2) ai z = L. In other words, 
the net variation over the extension A = fi2{L) — ^2{L/2) is A = —A. The net change of fi2 
over all the extended box is then zero. Thus the total charge of the extended system is also 
zero. The total charge however is a sum of the one inside the original box, Q, and that in 
the extension, Q. Thus Q = —Q. On the other hand, the charge Q can be calculated in the 
limit of slowly varying fj^2{z) from the formula p3|). Therefore one arrives at the equation 
(^) for the charge Q, independently of the details of the behavior of fi2{z) at intermediate z 
within the original bounding box. 



6 Charge density in sample models 



Here we consider explicit calculations of p{z) from the formula (^21) in sample models of the 
dependence of fi2{z) on z in order to illustrate that the distribution of charge generally does 
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not follow the slow variation limit described by Eq.(^) while the total charge is of course 
given by the universal formula (^. 

A limit that is maximally opposite to that of a slow variation is where /i2 chang es as a 
step function at a point in z. Here we consider the situation where [i2 changes at z = from 
a negative constant value /i2(-2) = —1^2 at z < to a positive one, ii2{z) = /i2 at z > 0. The 
Green's functions in Eg .(12^) can be found explicitly by the standard method of matching at 
2 = 0: 



1 



ptp 



1 

2^ 
1 

2^ 



1 + 



g-/i2 
q + fi'2 



exp 



^2q\z\) 
exp (—2 q \z\) 



(25) 



where q = yPo + /^i + /^i- The charge density is then given 'in quadratures' by the formula 



p(z) = /ii/i2 



B 
2^ 



oo 2tt 



-2^|pl + nl + nl \z\ 



(26) 



which, as one could expect, describes an exponential distribution around the discontinuity 
point z = over the range (/i^ + /i^)"^^^ that is the Compton wavelength of the fermion. 
For the total charge the integration over z and then over pq is elementary, and the result 
matches the general formula (^). 

It can be also noted, in connection with this example, that the non-exponential term 
in individual Green's functions in Eq.(p5D would lead to a divergence upon integration over 
Pq. However this term cancels in the difference, leaving the expression for p{z) finite. This 
illustrates the behavior discussed in general terms in the previous section. 

Another class of specific situations where the Green's functions are readily calculable in 
a closed form and the charge density can be found explicitly is presented by the Lee- Wick 



model In this model the term fi2 in the fermion mass arises from the interaction with 
a neutral field 0: P2 = 9 4>, with g being the Yukawa coupling. The 0^ self-interaction of the 
field 6 is assumed to lead to the domain wall solution 



(f){z) = V tanh(Az/2) 



(27) 



which interpolates between the two vacuum states (domains) with = ±v, and where A is 
the constant of the self-interaction, such that the mass of the bosons in either of the vacua 
is = Xv. 
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The behavior of the fermion states in such background depends on the ratio g/\ [|T5 
Clearly, the case g ^ X corresponds to the limit of slow variation of ^2, while the opposite 
case, g <^ X, corresponds to the approximation of an abrupt change in fi2, provided that in 
the latter case fii is also assumed to be small in comparison with m^. For a generic value of 
the ratio g/X the operators P^P and PP^ correspond to the solvable potentials of the form 
const/ cosh^(A z/2), and the Green's functions can be found in terms of hypergeometric 
functions. In a situation where 2 g/X is an integer: 2 g/X = N, both potentials are non 
reflecting, and the algebra is greatly simplified, since the relevant hypergeometric functions 
collapse to (Jacobi) polynomials in tanh(A2;/2) of the power N for P'^P and — 1 for 
PPK It can be also noted, that the case of = 2 in the considered class of models would 
correspond to a supersymmetric model, albeit with /ii = 0, which would take us beyond the 
assumptions adopted in the present paper. Here for illustrative purposes we consider only 
the most algebraically simple case of = 1. 

For X = 2g, the operator PP^ corresponds to a constant potential: PP^ = —d1 + (?^f^, 
and the corresponding Green's function is especially simple: 

exp(-g|2;i -Zsl) 



Z2 



2q 



(28) 



pl + fil + PP^ 

with q = \Jpq + fii + g"^ v^. The Green's function for the operator P^P can then be found 
16| , using the relations ([T0|) and (|l^) for its eigenfunction^: 

Uk{Zi) ul{z2) 



G(^PtP){zi,Z2; Po + ^1) = 
Vkizi) vl{z2) 



Zl 



1 



k=l 

1 



p.. 



pip 
Po + f^i 



Z2 



(29) 



P^ G^PPt){zi, Z2; 0) P - P^ G^ppi){zi, Z2; pI + 1^1) P 



+ 



Mo (^1)^0(^2) 



pI + fil 



81), and 



Applying in this formula the operators P^ and P as indicated to the expression 
using the zero mode from Eg. ([T^) : Uo{z) = ^gv/2/ cosh.{gvz), one finds that at coinciding 
points the Green's function ( pOj) is given by 

1 



where, again, q 



G(ptP){z,z; Po + /^i) 



Po + + 9^v^- 



+ 



2q 2q{pQ + fil) cosh {g v z) 



(30) 



^Thc individual Green's functions are finite, thus each can be expanded in the basis of the corresponding 
eigenfunctions. 
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Upon substitution of the results from the equations ( ]28|) and (^) to Eq.(p2D for the 



charge density, and after the integration over pq, the final expression for the distribution of 
charge in this model is found in a remarkably simple explicit form: 



where Q = B S aj:cta,n[gv / {f^l + g^v^)]/{27['^) is the total induced charge, in agreement with 
the general formula @. 

7 Discussion 

The explicit calculation described in the sections 2 and 3 illustrates the implementation 
of the general statement about the absence of spontaneous magnetization based on the 
C parity. In the concrete calculation the expected cancellation occurs due the fact that 
it is possible to regularize the theory while manifestly preserving the charge conjugation 
symmetry. Alternatively, as already noted, this cancellation can be viewed as the manifest 
Lorentz invariance of the regularized theory, since a nonzero spontaneous magnetization 
would certainly imply a breaking of the Lorentz symmetry In other words, neither 
the C symmetry nor the Lorentz one are anomalous. The same symmetries preclude an 
appearance of a spontaneous magnetization also in models where the fermions are assumed 
to have anomalous magnetic or /and electric dipole moments, including the cases [Q where 
the fermions are neutron-like, i.e. with zero charge and nonzero dipole moments. 

The really existing effect, allowed by both these symmetries, is the fermion charge density 
induced by an external magnetic field applied across the domain wall, which can be traced 
to the well-known anomaly in the axial current. The total induced charge is of a topological 
nature, and is determined, according to Eq.@, by the total magnetic flux through the wall 
and the asymptotic values of the fields in the model far from the wall. The distribution 
of the induced charge however, given by the general formula (0) depends on the details of 
the profile of the wall, as discussed in the sections 5 and 6. This effect, in principle, can 
be relevant in detailed analyses of phenomena at the walls in the early universe - either for 
flat domain walls, which could exist during some epoch, or for dynamics near the walls of 
bubbles during a first order phase transition. 

The arguments, based on the C parity and/or the Lorentz symmetry, do not directly 
apply to an asymmetric state with a net overall fermion charge, e.g. to the early universe 




(31) 
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with the baryon asymmetry, that is neither C nor Lorentz symmetric. In this case a sponta- 
neous magnetization of the wall is generally allowed, and is proportional to the asymmetry 
parameter of the considered state 0. It should however be clearly understood that even a 
magnetized domain wall does not produce a magnetic field, which is a simple consequence 
of the classical Maxwell's equations Therefore domain walls could not be a source of 
a correlated at cosmological distances magnetic field in the early universe, although their 
electromagnetic properties, like the charge distribution induced by a magnetic field, could 
generally be of importance in other aspects of dynamics of the early universe. 
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